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Theoretical analysis dealing with the interaction of electrons with the polar long-wavelength trans- 
verse optical (TO) vibrations is presented. The theory is based on the model of a polar crystal 
with classical potentials, which takes into account the electronic polarizability effects. A significant 
enhancement of the strength of the electron-TO-phonon interaction in ferroelectrics is found. A mi- 
croscopic justification of this effect is given. A bridge that relates the interaction of electrons with 
the polar long-wavelength TO modes of the lattice vibrations to the long-range dipole-dipole inter- 
action is established. As an application of our analysis, a new equation representing the relationship 
between the electron-TO-phonon interaction constant and material parameters is obtained. 
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O . I. INTRODUCTION 

C/3 . 

In the lattice dynamics of a polar crystal, there are two characteristic features, which give evidence that the 
interaction of electrons with the polar long-wavelength optical modes should be necessarily taken into account. The 
first feature is related to the existence of dipole moments associated with optical vibrational modes. The second 
one concerns with the difference in long-range fields given by longitudinal and transverse optical modes at the zone 
center In general, the long- wavelength optical vibrations in polar crystals are responsible for the presence of 
an electric field Q. As it is well-known, the generation of this field results from the displacement of the ions and is 
modified by electronic polarizabilities. This modifies the forces and affects the optical phonon frequencies. On the other 
hand, the anomalously large Born effective charges, which give giant LO-TO splittings in ferroelectric compounds 0], 
are associated with the existence of an anomalously large dipole-dipole interaction [3| . The long- wavelength optical 
phonons, therefore, appear to be a key factor for relating electronic and lattice (structural) properties. 

It is important to emphasize that electron-phonon coupling has the fundamental significance in the understanding 
£N| ■ of various physical properties of polar compounds (e.g., [j,[H|). However, to the best of our knowledge, most compre- 
hensive information regarding the interactions between electrons and optical phonons is primarily concerned with the 
theoretical and experimental studies related to longitudinal optical (LO) phonons (e.g., f4j u 9j and references therein). 
As regards transverse optical (TO) vibrations, with only a few exceptions [13, [Hj], no attempts have been made to 
consider consistently the problem of determining the electron-TO-phonon (el-TO-ph) interactions in polar crystals 
with due account for their particular polar properties. This question is of considerable practical importance for the 
vibronic theory (e.g., (l2rlla | and references therein) in which the coupling between electrons and the zone-center 
TO vibrations is the driving force of a ferroelectric instability, and where the strength of this coupling should be 
especially strong (l6| . Moreover, the new developments of the vibronic theory concerning different properties of oxide 
ferroelectrics and multiferroics [l3|, [TBI, IT7l - fl9j increased the interest in the understanding of the nature of the el- 
TO-ph interactions. Thus, this raises the issue of very little theoretical knowledge about the el-TO-ph interaction 
constants in polar compounds, and the present paper is aimed at filling this gap. 

In this paper, using characteristic parameters of a polar crystal in an explicit way, we study the main features 
of the el-TO-ph interaction. The novel aspects of our work are based on a first-principles treatment of the role of 
the el-TO-ph interaction. In particular, we demonstrate that, in polar materials, the el-TO-ph coupling is essentially 
influenced by the features of the Coulomb interaction between electrons and the lattice ions. We also give a microscopic 
justification of a significant enhancement of the el-TO-ph interaction strength in ferroelectrics compounds. Within 
a first-principles methodology, we show how to link the interaction of electrons with the polar long-wavelength TO 
phonons to the long-range dipole-dipole interaction. This, in turn, provides a nontrivial relationship between the 
el-TO-ph interaction constants and the macroscopic material parameters, such as dipole oscillator strengths and the 
forbidden gap. The relationship gives us a practical framework for estimations of the values of the el-TO-ph interaction 
constants in a polar material for which no data have been available so far. 
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In the more general context of understanding displacive structural instabilities our results can be used to probe 
deeper into the nature of the ferroelectric phase transition, and, in particular, to discover what factors related to 
material properties may be responsible for the strength of the el-TO-ph interaction in polar insulators. 



II. THEORETICAL MODELING 



A. Basic set-up of physical model and the main features 



As our focus is the el-TO-ph interaction, we shall be interested in finding the change in the electronic charge 
density induced by a polar TO vibration of a long wavelength. In order to consider the problem at the atomic level, 
we represent a polar crystal by a model based on classical potentials and contributions of the electronic polarization 
effects through the relative displacements of the electronic shells [l(| (see also 0, HJ). Within the standard procedure 
of the dipole approximation (expanding to first order in the ionic and electronic displacements, entering into the q- 
representation, and going to the normal coordinates u^j), the perturbation of the host crystal due to interactions with 
the polar lattice TO vibrations (qj) is given by: 



8U? xt {v)=N- l / 2 Y. V ^ v ) u <u « 



where Vqj(r) = uP(qj) • F q (r) , 
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P(qj) is an amplitude of the dipole polarization (dipole moment per unit volume) associated with the polar TO 
mode G denotes the reciprocal lattice vector, v is the volume of the unit cell, and N is their number. The 
derivation of Eq. (TTJ) accounts also for the fact that in the long-wavelength limit one can take the corresponding 
Fourier transformation of the Coulomb potential to be independent of ion's position in the unit cell. The accuracy of 
this approximation was investigated previously in [13, [H| • 

Equation (TTJ) is exactly the one that we shall use to study the interaction of electrons with the polar TO lattice 
vibrations in the long-wavelength limit. Before presenting the results themselves, let us emphasize the most important 
key points which characterize the model as a consistent, systematic and workable framework. 

The first is that Eq. ([1} accounts for the changes in ion properties that are caused by changes in its environment, and, 
hence, in accordance with [21| includes not only pairwise electron-ion interactions but also many-body effects. This 
follows from the accounting for an electronic contribution (P( e '(qj)) in the dipole polarization: P(qj) = (en) + 
p( e )(qj). This contribution is connected with the the distortion of the electronic charge distribution (e.g., [22|). 
The corresponding ionic contribution (PW (qj)) is related with the displacements of the ions participating in the TO 
vibrations [l[ : P^(qj) = v^ 1 ^2 s (Z s e/ \J M s )w(s, qj) where Z s e denotes the effective charge [23[ of an ion with the 
mass M s placed in the cite s, and w(s, qj) is the polarization vector of the normal TO vibration. At q = 0, the both 
polarization components satisfy the following relations 0, [l(| : 

P W (0j) = -^—P(0j) , P (e) (0j) = £ -^-^P(0j) (3) 

where eoo is the electronic (static-high-frequency) dielectric permittivity. 

Attention should next be directed to the long- wavelength nature of the quantity F q (r) standing in Eq. (fTJ). On 
the one hand, in a polar crystal it defines the electric field at a point r inside the bulk as a part of the internal field 
associated with the TO phonon displacement [Ic| (we recall that the microscopic electric field excited in the crystal 
has components q + G for all reciprocal-lattice vectors). On the other hand, at q = 0, F q (r)| q=Q = Fo(r) can be 
represented as (1 denotes lattice vectors): 

Fo(r) = -vf^^-') =-We(r). (4) 



It is seen from Eq. (0} that, in accordance with [24|, H^, the quantity Fo is a field generated by the Coulomb 
electrostatic potential e (r) which is connected (at given boundary conditions) through the Poisson equation with a 
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charge distribution 6p e (r) = ej^|£(r— 1). This result, in turn, suggests a macroscopic characterization for the quantity 
Vqj(r) given in Eq. ([TJ as follows: it conies in the form of an electrostatic contribution which interrelates the dipole 
polarization associated with the polar long-wavelength TO mode and the electric field caused by the reaction of the 
electron subsystem to the relevant "external" perturbation. In this connection, we also note that both the resulting 
shift of the electron density and the effect of the dipole polarization appear to be clear signatures of the strength of 
the el-TO-ph interaction. 

The third important point is that interaction in the form (JTJ) may be responsible for the softening of the lattice 
subsystem. Clearly, this is the case of the vibronic theory 12|, H| • Following closely Ref. [26| , we observe that local 



change in the electron density distribution will generate an imbalance with respect to the corresponding density of the 
positive background. As a result, there will be a restoring force that will tend to compensate this density imbalance. 
As it argued in (26|, such asymmetry effect will initiate distortions of the ions and the redistribution of charge, 
eventually leading to a decrease in the total energy, and to a match between the symmetry of the ionic arrangement 
and that of the electron density distribution. 



B. Model for electron-TO phonon interaction 



The starting point is a generic model of an insulator consisting of the occupied valence and empty conduction 
bands separated by a forbidden gap. Keeping in mind correctness in the description of an insulated state, we assume 
that the wave functions (ipak) and energies (E a (k)) are the set of one-particle Bloch eigenfunctions and eigenvalues 
corresponding to the band structure, which are determined within the framework of DFT-LDA methods with the 
correct accounting for the quasiparticle effects by means of many-body perturbation theory [27| . 

Equation (JTJ) , providing a basis for studying the el-TO-ph interaction, corresponds in the second-quantization 
formalism to the Hamiltonian 

H el _ ph = N- 1 ' 2 J2 J2 9aa' ( k > <U') <£kVk-q «qj (5) 
a, a j kq 

which characterizes the dynamic mixing of an electron state | a k > with another state \ a (k — q) > caused by the TO 
phonon. Here a + (a) are the creation (annihilation) operators for electronic states in the valence and conduction bands 
(a, a ), g aa ' (k, qj) denotes the matrix elements of the el-TO-ph interaction defined at the equilibrium high-symmetry 
configuration of the lattice by 

<w ( k > <u) = < ^\ v<u Wk' (k - q) > ■ (6) 

Expression (jB]) involves both interband (u 7^ a ) and intraband [a = cr ) matrix elements. In the long-wavelength 
limit, the intraband matrix elements vanish due to inversion symmetry and, correspondingly, the coupling of electrons 
and TO lattice vibrations is completely characterized by the interband matrix elements. 

In the context of first-principles features/attributes belonging to the model, it is worth to emphasize the following 
aspects: (i) The electronic states are periodic in space and delocalized over the entire unit cell; in systems with a 
gap infinitesimal displacements affect all the states uniformly |28[, (ii) the valence and conduction (excited) states 
can mix under these displacements [II], (iii) the energetics of the interaction of electrons with the long- wavelength 
polar TO lattice vibrations is strongly affected by energetics of the lattice electronic distribution, (iv) the effect of the 
dipole polarization is to change the center of force between a "shell" of an ion's electrons and the other force centers 
in the crystal [29] , (v) like the electron-gas model of Gordon and Kim [3(| , the corresponding interaction energy is a 
function of changes in the electron charge densities, and (vi) as seen from Eq. these changes are represented by 
distortions of charge density owing to electric field gradients at certain lattice sites of a polar crystal. 

Our intermediate goal here has been to offer a detailed, microscopically realistic framework that could be used for 
quantitative investigations of the el-TO-ph interaction in a polar crystal. In order to achieve the objectives of the 
present work further study in this direction is divided into two main steps: The first is a comparative first-principles 
analysis of the contribution of Eq. ([5]) into the lattice dynamics of TO lattice vibrations; the second handles the 
interband matrix elements, and develops formulas allowing us to estimate el-TO-ph interaction strengths in terms of 
macroscopic material parameters for a wide range of polar dielectrics. 
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III. ELECTRON-PHONON INTERACTION AND TO LATTICE VIBRATIONS 
A. On the dynamics of TO lattice vibrations 

We begin the present section by reviewing the main properties of the electronic contribution into the dynamics of 
lattice vibrations. Then we calculate the relevant electronic contribution given by Eq . ([51) . Our task is to compare 
this contribution with the corresponding results of microscopic lattice dynamics |3ll433j with regard to the long- 
wavelength polar TO phonons. For this purpose, within the adiabatic approximation, we consider the dynamics of 
the TO vibrations in terms of the matrix elements of the el-TO-ph interaction determined by Eq. ([6]). 

A key feature of the electron subsystem is that it provides a link between an electronic energy and an equilibrium 
geometry of a crystalline structure (this link in turn depends upon the electronic contributions to the total energy). 
The electronic contributions to the dynamical matrix of a nonmagnetic crystal arise from changes in the electronic 
charge density due to the presence of the electron-ion potential in a system (e.g., [HI, 0, [3l]> [34], HH). This implies 
that, by studying the lattice dynamics, one can acquire fundamental information on the role of the electron-phonon 
interactions in a given material. 

Since the real part of the self-energy of the TO vibrational mode specifies the deviation of the unperturbed phonon 
states upon adiabatic variation of the potential connected with external perturbation, we represent the square of the 
TO phonon frequency Q 2 ^ as a sum of the unperturbed term and the additional term Aw^- which accounts for 
cl-TO-ph coupling governed by Hamiltonian (J5j: 

We expect the "normal" part w^, which may be regarded as the bare TO phonon frequency with respect to the 
el-TO-ph coupling, to be free of softening anomalies. As usual (e.g., [HI, HU), it is assumed that uj 2 ^ also involves 
contributions from the phonon-phonon interactions, which are needed to describe the dependence of the TO phonon 
frequency on temperature. 

We can obtain the corresponding electronic contribution Au^- in the following way. The first step is to employ 
the linear density response function matrix (the dielectric susceptibility) \ that relates a variation of the electronic 
density (the charge response) Sp to the perturbation potential Sv ex t as follows [37l[38j: 

8p = x5v ext , e~ 1 =I + v c x (7) 

where e _1 denotes the inverse dielectric matrix, / is the identity matrix, and v c is the bare Coulomb interaction. Note 
that perturbation (J5J) modifies the ground-state values and therefore induces Sp. As stated by the DFT formalism, 
both x & n d the potential Sv ext itself are unique functionals of the electronic density. 

Since the Coulomb and electron-phonon interactions cause the phonon propagator to be fully renormalizcd 
(e.g., 39]), the second step is to represent the real part of the TO phonon self-energy to the lowest-order in g 2 
as the corresponding harmonic term, and to account for the contributions of the Coulomb interaction among elec- 
trons. Using microscopic expressions for matrix elements g acr '{k., qj) given by Eq. ©, by the last step, we write the 
equation for Au^- in the following form: 

= 4nv B a p (q) P a (qj) Pp (qj) (8) 



where 



r 



B Q/3 (q) = ^ ]T G a v c (q+G) Xe i(q + G,q + G)v c (q+G)G 0) (9) 
g,gVo 

** = r ^fT, vM + G) = ^— 2 , (10) 
I-v c Il v|q + G| 



{ ' iv^4^(k+ q )-£ ff (k) 



<7.tT 



J C k e l(q+G)r VVk- q ^ / C' k - q ^ (q+GV '^dr' • (11) 



Here the matrix elements of the microscopical dielectric susceptibility of the electron subsystem Xel{q+ G,q + G ) 
are expressed in the random phase approximation in terms of polarization operator II (37l . l40l - |4^ | , and / are the 
occupation numbers. 
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B. Consistency with first-principles considerations 

In this and next section, we shall show that our results given by Eqs. (|8j)- (fTTj) are consistent with first-principles 
considerations. This will allow us to establish a bridge that relates the interaction of electrons with the polar long- 
wavelength TO modes of the lattice vibrations to the long-range dipole-dipole interaction. 

Within the framework of our comparative analysis, it seems to be useful to take into account three important 
points: 

(i) We deal with the long- wavelength TO phonons, i.e., in fact, with the analytical part of the harmonic force 
constants. The remaining non-analytic part only affects the LO phonon frequencies, and therefore determines the 
LO-TO splitting at the T point. 

(if) The unique feature of the force constants is that they can be decomposed into a sum of two independent terms [3l], 
HjjUE]: the direct ionic and electronic contributions. The electronic contribution, which represents indirect interactions 
via the electron subsystem, is the focus of our discussion because involves the detailed information on the el-TO-ph 
interaction. 

(iii) The long-wavelength behavior of the force constants can also be represented in terms of the balance between 
the short-range repulsive and long-range Coulomb forces associated with the short-range and long-range dipole-dipole 
interactions, respectively [1, |32|. l33l |43 |. The long-range contribution is represented by the interplay of the Born 
transverse effective charge tensors. It is a key quantity in polar compounds because favors the ferroelectric distortion, 
while its direct competitor (the short-range contribution) tends to suppress the ionic displacements and to provide 
the stability of the high-symmetry configuration. 

Note that, in the context of the present work, the above pair- wise separabilities of the force constants will be 
particularly instructive for comparison purposes. For instance, by matching different pieces of information contained 
in the dynamical matrix, we could recognize relevant aspects of the el-TO-ph interaction. 

We start our comparison procedure by noting that within the microscopic theory of lattice vibrations the Fourier 
transforms of the harmonic force constants are expressed through the matrix elements of the susceptibility function 
Xei(q + G, q + G ) (e.g., [3l|, H3, El])- By performing a comparison of Eqs. © and © with those from [3l|, we can 
ensure that the expression for _B Q( g(q), which is written in terms of the microscopic electronic dielectric susceptibilities 
and the Coulomb potentials, describes the corresponding electronic contribution to the square of the TO phonon 
frequency fi^- (i.e., the contribution which uniquely involves the Xei(q + G, q + G ) matrix elements with both G 
and G different from zero). Due to the microscopic formulation of the el-TO-ph interaction given by ©, this result 
is general and independent on the particular details of the interaction of electrons with the polar long- wavelength TO 
vibrations. This implies that fundamental property of B a p(q) is that it represents, in accordance with [3ll l39l. |45|. 
the electron mediated part of the long-range Coulomb interaction. 



C. Relevance to the long-range dipole-dipole interaction 

As the second step of the comparison procedure, it is a ppr opriate to compare equations (JH) and ([9]) with the 
corresponding results of lattice dynamics theories [l], H^, I3ll-l33j related to polar insulators. For this purpose let us 
inspect the long-wavelength limit of Eqs. ([5]) and © in more detail. 

Observe that one of the characteristic features of Eq. (JSJ represents the interplay of the polarization amplitudes 
Pa(0j) and Pp(0j) which, as in [33j, can be characterized by the second rank tensor S a p(j) = 4ttv P a (0j)Pp{Qj). The 
quantities S a p(j) are defined as dipole oscillator strengths for the given zone-centre TO vibration and correspond to 
the contributions of the TO vibrational modes in the infrared (IR) part of the Lorentz model of the dielectric matrix 
e a p (e.g., EH 13): 



(e w ) Q /3 = (£oo)q/3 + ^2 r,2 aP _ 2 • ( 12 ) 

s L 0j w 



It is easy to show, using the definition of Pa\0j) together with that the dipole oscillator strengths S a p{j) can 
be represented in terms of the TO mode effective charge: 

SapU)=^- \ >Z. i I > Z, — =^ \ . ,!:>>) 



Here the result of each such summation in the round brackets is called the TO mode effective charge |2| (which 
corresponds to the dipolar activity of the zone-center TO phonons and is a measure of the intensities of IR-active 
modes), Z* are the Born transverse dynamical effective charges defined as [l], HH, H3, H§| Z* = {(too + 2)/3] Z s . 
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After substituting for the product P a (0j)Pp(0j), Eq. <j8j) becomes 



Aa4 = J2 B <xp(Q)S<*p(j) (14) 



where the tensor B a p(Q) is defined from Eq. © at q = by 

B «rf(°) = 4^2 E G a « c (G) Xei (G,G> c (G')G^ (15) 
g.gVo 

and x e i(G, G ) = lim q _i.o Xei(q + G, q + G ). By using the representation of the susceptibility function Xei(G, G ) in 
terms of matrix elements of the microscopic polarizability tensor a 7?? (G, G ) (49[, 

Xei(G, G ) = - ^2 G 1 a 7n (G, G ) G q , 

Eq. (fl~5|) becomes 

*M°)=-^ E E^^ G ')§i- ("J 

The set of equations (|14[) - (|16p . which we have rigorously derived, allows us to analyse the relevance of our model to 
long-range dipole forces. We shall now prove this relevance by comparing with the results of the previous theoretical 
studies of lattice dynamical models (e.g., [3, US ISll)- 

Note first that, in polar insulators, in the limit of zero wave-vector the displacement of charges from their equi- 
librium positions creates dipoles, which interact with long-range forces [l|. The associated long-range dipole-dipole 
interaction can be characterized by the coupling of the Born transverse dynamical effective charges Z* Z£, i.e., by 
the product P a (0j) Pp(0j) . Secondly, this interaction can be characterized by assigning the dipole polarization an 
external moment P(0j) = [(coo + 2)/3l PW(Oj') with the proportionality coefficient given by the local-field factor 
(see Eqs. (J3|)). Thirdly, following Ref. [5lJ, Eq. (TT4"]) with the help of Eqs. (fT5|) and ([T5]) can at once be identified as the 
long-range Coulomb contribution resulting from the induced dipole interactions through the electronic polarizability. 

Therefore, our comparative analysis shows that the set of equations (fT4| -(fT6 |) are directly related to the contribution 
of the long-range Coulomb interaction. 

Furthermore, from comparison of our result given by Eq. (|14p with the corresponding equations for the TO vibra- 
tional modes derived within the macroscopical formalism we can find the evaluation of the tensor B a p(Q). Let us 
show that B a f; (0) is expressed as 

B ol0 {Q) = -8 ol0 /{e oo +2). (17) 

As a comparison parameter, it is reasonable to choose the product P a (0j) Pp(0j) that measures the long-range dipole- 
dipole interaction [Ho| • With the use of results of 0, , for the case of the cubic (or tetrahedral) lattice symmetry, 
in dipole approximations, the counterpart of Eq. (fM)) can be represented in terms of the product P a {0j) P/3(0j) (or 
Sap{j)) as follows: 

Awgj = -4nvJ2-^P«(0j)Pe(0j) = -E;%^(J)- (18) 

Eq. P7|) is obtained by a direct comparison of Eqs. (TT4"|) and (TT8"|) . 

In the last step of our comparison analysis, we find that equations (fT4|) and (fl7|) are identical to the results derived 
within the exact accounting of the dipole-dipole interaction in lattice dynamics [33j . This identity directly implies 
that the product B a p(Q)P a (0j)Pf)(Qj) represents the regular at the q = contribution of the long-range dipole forces. 
This allows us to conclude that, within the framework of TO vibration mode dynamics, the model originally specified 
by Eqs. <|T]) , ((2j) , d5j) and (j6|) matches the microscopic approach proposed and developed for polar crystals in [32l.l33|. 
This is a novel result of principal importance: we showed and confirmed the correspondence between two descriptions 
of dynamics of the zone-centre TO vibrational mode - one in terms of the el- TO-p h interaction and the other in 
terms of the long-range regular at the q = part of the dipole-dipole interaction [32l |33||. Moreover, Eq. (fT4"|) clearly 
demonstrates the equivalence of the interpreting power of both descriptions, and therefore can serve as a link between 
the two approaches. 
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Thus, we have shown above that microscopic formulation of the el-TO-ph interaction given by ([6]) provides results 
consistent with first-principles considerations. We finish this section with a few important remarks. Firstly, it should 
be noted that, as follows from the foregoing analysis, the interaction of electrons with the polar zone-centre TO 
phonons is directly associated with the long-range Coulomb interaction. Secondly, the magnitude of the destabilizing 
contribution (B a p(0) < 0) determined by Eq. (fl4l) will rise steeply in polar materials. This will lead to significant 
softening of the zone-centre TO vibrational mode frequency J7oj ■ The third remark is concerned with strong sensitivity 
of calculated (at q = 0) frequencies ([TJ| to the dipole oscillator strengths: this is the manifestation of the fact that 
the zone-centre TO vibrational modes exceptionally strongly can be coupled to valence band electrons. 

Our findings give a clear understanding of the role of the el-TO-ph interaction in the conversion of f2(y to the 
ferroelectric soft-mode. Eq. ([Til) shows how changes in the electron density (or the changes in the electronic states) 
induced by the local displacements are sufficient to yield a considerable softening of the TO phonons around the 
r point. In the context of structural dynamical instability the dipole oscillator strengths S a f3{j) as enhancement 
factors in Eq. (|14l) can serve as an important macroscopic measure for probing how close the system may be to the 
ferroelectric instability. 

In the next section, we specify our treatment of the el-TO-ph interaction to the case of ferroelectric crystals, 
and derive the corresponding effective Hamiltonian, which would characterize the interband scatterings of the bond 
electrons due to the zone-center TO phonons. 

IV. ON ELECTRON-TO-PHONON COUPLING IN FERROELECTRICS 
A. Parametrization of electron-TO phonon coupling at q = 

The microscopic approach considered above is quite general; in the present section, we show how simplifications 
allow us to deduce a reduced model of el-TO-ph interaction. 

The starting point is Eqs. (|8|)- dTT1) in which we set q = 0. The main complication in executing the limit q — > 
in the matrix elements of \ei is that one needs to exclude the relevant contributions associated with the electronic 
displacements (these contributions are included in the definition of g 2 ,). This problem is solved in two steps: (i) 
by finding a change of the electronic charge density associated with interaction ((TJ) and (ii) by applying Eqs. ([7]) in a 
similar way as was described in [32] regarding the analysis on the internal electric field effects. Acting in this manner, 
and taking into account Eqs. ([3]), we obtain that in the long- wavelength limit, the reduced susceptibility \ei turns out 
to be related to Xei by the following equation: 

q -> : Xei = — ti ■ (19) 

Coo T" £ 

With all of these considerations and assuming that in insulators (/ CT (k) — f a '(k)) = 1, we can write the long- 
wavelength limit of Eqs. © in terms of the el-TO-ph interaction g aa > (k, 0j)=g 3 , (k) as follows: 

Ao, 2 - 3 1 V V l^'W 2 (on) 
AuJ V--7^T2N±; \ |^(k)-^(k)| - (20) 



l^(k)| 2 = ^E P «(°i)^(°j) E G a v c {G)v c {G)G p 
a, p g.g'^o 

x J C k e lGr Wk^ / VC< k e- iGV V CT kdr. (21) 

The remarkable property of Eq. (|2TJ|) is that it can be represented as a product of two independent factors, (A- 7 ) 2 x 

flo(0), one of which, (A-?) 2 = | g 3 | 2 [{too + 2)/3] _1 , characterizes an effective interaction of electrons with polar zone- 
centre TO vibrations. The other, n (0), is the so-called non-interacting susceptibility. Following the vibronic theory 
prescription (e.g., [121] 1. this allows us to constitute a reduced one-parameter version of the el-TO-ph Hamiltonian ([5]) 
in the following model form: 

= £ £ A^ (k) a+ k(Vk u 0j . (22) 
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First-principles justification of Eq. (|22|) is an important result of the present work. We refer to this equation as 
a vibronic-type model since it involves the interband coupling between electrons and the zone-centre TO phonons. 
According to [f2|, El] , such coupling forms the basis of the vibronic theory. Due to the successful mapping of the 
original microscopic model ([5]) onto the effective model given by (|22p . it became clear why, by introducing the constants 
X J as model parameters, it is, in principle, possible within the framework of the vibronic theory to study the lattice and 
electronic properties of polar materials. Besides, Eq. (|20[) shows the new meaning of the el-TO-ph coupling constant, 
the square of which is inversely related to the local field factor (e^ + 2)/3 . Since local fields and their effects become 
more prominent in compounds with a mixed ionic-covalent character of chemical bonds (e.g., [371. |52|. |53||). this would 
certainly imply the important improvement in the applicability of the standard vibronic Hamiltonian. Note that a 
similar renormalization occurs in the theory of surface polarization modes || . 

From a physical point of view, Eq. (|22[) represents the Frohlich-type long-range electron-lattice dynamic hybridiza- 
tion of the electronic bands of opposite parities, which takes into account the relevant s-, p- and ci-channels (e.g., [T3 - fl5| 
and references therein). For instance, as applied to a family of high polar crystals such as the ferroelectric ABO3 
perovskite oxides, hybridization (|22|) involves significant mixing between the filled O 2p and the em pty d° (Ti 4+ , 
Nb 5+ , Zr 4+ , Ta 5+ , Mo 6+ , W 6+ , etc.) electronic states caused by the IR-active TO Fi„ soft vibrations flS. HI. l53| . 

For the last several decades, the model of t ype (122p was a useful prototype for intensive studies of various properties 
of both typical insulating perovskites [l3l . \vK 15414571 ] and the A IV B VI narrow-gap semiconductors and their alloys [58l — 
[62J. It was recently demonstrated to be applicable to give a theoretical description of the ferroelectricity found in 
BiFe03-type multiferroics [18|, Hlf ■ 



B. Evaluation of electron-TO-phonon coupling at q — 



In order to evaluate the el-TO-ph coupling in the long wavelength limit, we first examine the Frohlich-type expression 
for the amplitude of the interband el-TO-ph interaction g atJ i (k, qj). This can be accomplished in a similar way as 
in the case of the scattering of an electron due to LO phonons [lOj. Let us fix j. Denoting the amplitude by 
I 5 CT(T '(k)| = (liniq-^o I g aa ' (k, qj)! 2 ) 1 / 2 and using the evaluation of the dipole oscillator strengths of the TO modes 
Sapij), we can obtain: 



|<w(k)| - An LT 



4-7T 



< crk|F (r)|cr'k > 



(23) 



where, according to [63j j . we introduced for the given polarization a difference between the longitudinal (Qlo) and the 
transverse (flro) frequencies (the LO-TO splitting) provided in terms of ASI^t = {£l 2 LO — fif^) 1 / 2 (i.e., the difference 
in long-range fields given by longitudinal and transverse modes as q — > 0). 

From Eq. (|23j). one can draw the following important inferences: First of all, note that dynamic hybridization 
described by ([5]) leads to an asymmetric charge distribution, which corresponds to the internal electric field Fo(r). 
This field, in turn, provides, due to Eq. ((4]), the long-range character of the el-TO-ph coupling at the T point. Secondly, 
the magnitude of the splitting Ail lt serves as an enhancement factor of the matrix elements of Fo (r) . Therefore, in 
the long-wavelength limit, the el-TO-ph interaction in polar crystals is distinguished from that of most other dielectrics 
by the following features: (i) it is long-range, (ii) it is controlled by the internal electric field and (iii) it is essentially 
sensitive to values of the LO-TO difference (i.e., the polar strength of long- wavelength optical vibrational modes). 
Eq. (|23l) can serve as indicative of the strength of the el-TO-ph coupling in polar materials. In particular, large values 
of I g aa ' (k)| occur when AH^t is large: for example, the experimental difference between the corresponding LO and 
TO modes in the ferroelectric BaTi03 is about 530 cm -1 [64[ to be compared with the values of 100, 69, and 48 cm -1 
obtained for the differences in the compounds NaCl, KC1, and RbCl, respectively (65[. 

In a polar lattice, the LO-TO splitting depends generally on Born's transverse dynamical effective charge Z* of the 
lattice ions and the screening of the Coulomb interaction, which depends on the electronic dielectric permittivity £00 Q : 
An 2 LT cx |Z*| 2 /£oo- Thus, one can conclude that the presence of anomalously large Born's effective charges is the key 
signature that the interband el-TO-ph couplin g is essentially strong in a ferroelectric material. This result corresponds 
to the main assumption of the vibronic theory [12| - |l5| regarding the existence of sufficiently strong interband el-TO-ph 
coupling in displacive ferroelectrics. The strength of the interband el-TO-ph interaction can therefore serve as a direct 
indicator of the extent to which a crystal lattice is close to a possible ferroelectric instability. Further details concerned 
with the quantitative characteristics of the electron-TO-phonon coupling are given in our separate paper [l6| . 
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C. Electron-TO-phonon coupling via macroscopic parameters 

As shown in the present work, one can associate the interaction between electrons and the polar long- wavelength TO 
phonons with the long-range dipole-dipole interaction. This link, which was not anticipated by previous theoretical 
considerations, helps us to express the strength of the el-TO-ph interaction via macroscopic parameters of a polar 
crystal. 

Let us introduce the bare constant g j which represents the effective k-independent interband el-TO-ph interaction 
defined by a proper average of the squared el-TO-ph matrix elements over the given electronic states: 

<=Y.^Y.Kjn \ EA ^'_ EAn (24) 

where E g is a bond-gap energy. Accounting for local-fields effects-induced partial screening of the bare el-TO-ph 
interaction is described by the renormalization: 

9oj — > Sj = 3^/(600 + 2) . (25) 

As follows from Eqs. ([T4|) and (fT?]) . the contribution of the linear el-TO-ph interaction to the square of the zone- 
centre TO vibrational mode frequency f2g ■ can be represented as the product of two macroscopic quantities: the 
quantity B a p(0) which accounts for the corresponding electronic contribution, and the dipole oscillator strengths 
S a f}{j) associated with the given TO vibration. The bare constant of the el-TO-ph interaction goj can thus be 
determined from matching Eq. ([20)) and Eq. (fTl| together. As a result, we obtain the following relationship (S(j) = 

£ a SaaO')) : 

4 • = ^E g MjS{j). (26) 

Eq. is an important result of the present work. It relates, for each zone-center TO vibration of the branch j 
and of the reduced mass Mj, the bare constant of the el-TO-ph coupling at the T point, goj, with the macroscopic 
material constants, the values of which can be obtained from experiment. Employing the relevant information from 
the IR spectra concerning the IR-active optical phonons and the dielectric function behavior in the far-IR spectral 
range, and using Eqs. (I26|) and (I25[) in combination with the experimental data for and E g , we can directly evaluate 
the el-TO-ph coupling constants for polar materials of interest. 

A more detailed demonstration of the practical usefulness of Eq. (1261) , together with numerical results for a wide 
number of selected polar insulators and semiconductors, can be found in our separate publication jl6|. In particular, 
calculations of the interband el-TO-ph interaction constants and the further comparative analysis showed that the 
large interband el-TO-ph interaction is a special microscopic feature of the ferroelectric materials. In contrast, in 
non-ferroelectrics, as it was demonstrated, the strength of the el-TO-ph interaction is not necessarily high enough due 
to their lower polar nature. 

V. CONCLUSIONS 

In this paper, focusing mostly on the fundamental contribution of electron subsystem to the dynamics of polar 
long-wavelength TO vibrations and using a first-principles methodology, we provided a systematic description of 
the el-TO-ph interaction in a polar insulator. The study is based on the model of a polar crystal with classical 
potentials, which takes into account the electronic polarizability effects. By analyzing the electronic contribution to 
the TO vibrational mode in terms of the el-TO-ph coupling, we established a bridge which allowed us a) to link the 
model under consideration to the microscopic lattice dynamics and b) in the long wavelength limit, to relate the 
interaction of electrons with polar TO vibrations to the long-range dipole-dipole interaction. Our results highlight 
the importance of the el-TO-ph interaction for the genesis of the long-wavelength TO vibrations in ferroelectrics, 
thereby giving fundamental support at the microscopic level for the applicability of the vibronic theory. Within a 
first-principles methodology, we found and explained the significant increase of the constants of the el-TO-ph coupling 
in ferroelectric materials, showed how the el-TO-ph interaction constants might be dependent on macroscopic material 
parameters, and obtained analytical equations allowing us to estimate el-TO-ph interaction strengths in a wide range 
of polar dielectrics. In particular, it was proved that the zone-center TO vibrational mode effective charge Z* can be 
considered to be the key macroscopic parameter of the el-TO-ph coupling strength. In materials where the el-TO-ph 
coupling is operative, it can be verified by spectroscopy measurements of the IR-active TO mode. 
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